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Abstract 


:£^he planetary exospheres are poorly known in their outer parts, since the neutral densities are low com- 

O 

Chared with the instruments detection capabilities. The exospheric models are thus often the main source of 
^Iformation at such high altitudes. We present a new way to take into account analytically the additional 
Effect of the stellar radiation pressure on planetary exospheres. In a series of papers, we present with an 




Hamiltonian approach the effect of the radiation pressure on dynam ical trajectories, density prohles and 




scaping thermal flux. Our work is a generalization of the study by 


Bishop and Chamber laini fjl989l) . In 


^his third paper, we inyestigate the effect of the stellar radiation pressure on the Circular Restricted Three 

O, 

^ody Problem (CR3BP), called also the photograyitational CR3BP, and its implication on the escape and 
’^he stability of planetary exospheres, especially for Hot Jupiters. In particular, we describe the transfor- 
Ljiaation of the equipotentials and the location of the Lagrange points, and we proyide a modihed equation 
^^r the Hill sphere radius that includes the influence of the radiation pressure. Finally, an application 
^o the hot Jupiter HD 209458b reyeals the existence of a blow-off escape regime induced by the stellar 
Radiation pressure. 

'^^k.eywords: exosphere, radiation pressure, CR3BP, escaping flux 

OO 

o- 

’—1. Introduction 


• 

^ The exosphere is the upper layer of any planetary atmosphere: it is a quasi-collisionless medium where 

^5-h 

cUhe particle trajectories are more dominated by gra vity than by collisions. Ab oye the exobase, the lower 
limit of the exosphere, the Knudsen number Kn (IFerziger and Kaperl. Il972h becomes large, collisions 
become scarce, the distribution function cannot be considered as maxwellian anymore and, gradually, the 
trajectories of particles are essentially determined by the gravitation and radiation pressure by the Sun. 


The trajectories of particles, subject to the gravitational for ce, are completely solved with 


of motion, but it is not the case with the radiation pressure (IBishop and Chamberlainl. Il989h . 


he equations 
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^ pressure disturbs the conics (ellipses or hyperbolas) described by the particles unde? 

the influence of gravity. The resonant sc attering of so l ar ph otons leads to a total momentum transfer 
from the photon to the atom or molecule (IBurns et ahl. 119791) . In the non-relativistic case, assuming an 
isotropic reemission of the solar photon, this one is absorbed in the Sun direction and scattered with the 
same probability in all directions. For a sufficient flux of photons in the absorption wavelength range, the 
reemission in average does not induce any momentum transfer from the atom/molecule to the photon. The 
momentum variation, each second, between before and after the scattering imparts a force, the radiation 
pressure. 


Bishop and Chamberlain (119891) analyzed its effect on the structure of planetary exospheres. Neverthe¬ 


less, their work was limited only to the Sun-planet axis, with a null component assumed for the angular 


momentum around the Sun-planet axis. We thus generalize here their work to a full 3D calcu 


order to investigate th e influence o : 


as the density profiles flBeth et al. 


the ra diation pressure on the trajectories flBeth et al. 


ation, in 


2ni5al) . as well 


2015bl) and escape flux (future work). 


In this paper, we propose to investigate the influence of the radiation pressure on the stability of 


planetary exospheres, based on dynamical effects of t 


validity of our approach (see. 


Beth et al. 


le radiation pressure. In particular, we study the 


(l2015al. l2015bl) for the specific case of Hot Jupiters, which needs 


to consider the Circular Restricted Three Body Problem (CR3BP) combined with the radiation pressure 
i.e. the photogravitational Three Body Problem. We will also derive a new expression of the Hill’s sphere 
radius (the limit of the gravitationnal influence of the central body) that includes the effect of the radiation 
pressure, before we discuss the impact of the radiation pressure force on the stability of the Hot Jupiters 
atmospheres (with a case study on HD 209458b). These may be indeed be bounded or not to the planet 
depending on the radiation pressure intensity. 

In the section [21 we develop our approach to study the photogravitational CR3BP. Then, in the section 
131 we will present the topology of equipotentials, derive a new expression of the Hill’s sphere and then 
discuss the possible consequence for the planetary atmosphere stability and the escaping flux. Finally, the 
conclusions will be reported in the section 01 


2. Formalism 

We investigate here the combined effects of the external forces on planetary neutral exospheres: the 
planetary gravity, the stellar gravity, the centrifugal force and the radiation pressure. As we know, the 
radiation pressure is directly proportional to the flux of photons. Neglecting relativistic effects such as the 
Doppler shift or the Poynting-Robertson effect, the radiation pressure then depends only on the square 
of the distance from the star as the stellar gravity. We define the 13 parameter as the ratio between the 
acceleration induced by the radiation pressure and by the stellar gravity, this parameter being a constant 




































tEP^MA^^f^^ystem. Indeed, the radiation pressure effect without relativistic effects depends only o3 
the flux and thus the square of the stellar distance as the stellar gravity. 

Each neutral is then subject to the dimensionless effective potential —r2(a:, y, z) in the rotating coordi¬ 
nate system defined as: 


(-)/ \ ^ / 2 , 2'i I P) I 

= -[x +y) + -+ — 


rs 


Tpl 


( 1 ) 


and 


rs = -s/^x + yY + y‘^ + z"^ 

Tpi = ^/{x + y-iy+ y^ + z"^ 

with y = Ml/{Ml -|- M 2 ), Mi the mass of the planet, M 2 the mass of the star, (—/i, 0 , 0 ) the position of 
the star and (1 — p, 0, 0) the one of the planet. The axis are oriented such as the x-axis is the star-planet 
axis, the z-axis is perpendicular to the planetary motion and the u-axis so that t 
orthonormal basis. The notations are similar to the ones used in 


y-axis so that t 

le Oxvz frame is a direct 

Simmons et al. 

(19^5 

)• 


The first term of hi describes the potential energy associated to the centrifugal force, the second one 
corresponds to the stellar gravity and the radiation pressure effect 13 and the last one the planetary gravity. 


by 


This problem was prev iousl y studied concerning 


Schuerman (1972 


1980) and 


Simmons et ah 


h e positions and th e stability of equilibrium points 


(119851) . ISchuermanI (jl972h determined the position of the 


Lagrange points L 4 , L 5 , then investigated the equipotentials and the (3 values for which the Roche lobe of 
the planet is connected to the ste’ 
both are not connected at all (cf. 


Schuerman 

(I 972 

ao a iuiii_.Liuii ui /x cliiu. /j. ± k 

), figure 4). On another hand. 

JL VCl V DliiCLll 0/ 

Simmons et al. 

(19^5) 


investigated the general case of a binary star system, where each star has its own radiation pressure. They 
gave a complete review about the Lagrange points Li, L 2 and L 3 with the numeri c al der ivation of their 
positions and their stability, whereas L 4 and L 5 were mostly studied by 


Schuerman] (119801) . 


The link between t 


et ah (201 


le stability of a planetary atmosphere and the CR,3BP was tackled bv Lecavelier 


des Etansrs et al. (120041! about the strong escape rat e observed on HD 209458b and wide l y cove red in the 


liter ature (IVidal-Madiar and Lecavelier des Etangsl (12004) 


Lecavelier des Etangs et al. 


[ 2004 ^Koskinen 


)). They modelled with different densities and temperatures estimations the upper atmosphere 


of HD 209458b. For high thermospheric temperatures, the thermosphere could extend so far that the 
exobase could reach the Roche lobe limit, located at only ~ 4 planetary radii. Thus, in this case the 
requisite kinetic energy to escape from the exobase (or critical radius r^) is smaller than the common 
gravitational potential GMi/tc- Several papers dedicated to HD 209458b investigated this issue, focusing 
on the mass loss rate and the conversion of the EUV flux into kinetic energy for Hydrogen so that it 
can escape. However, most of the previous studies did not take into account the additional effect of the 
radiation pressure (i.e. /3 = 0 ) which acts as a repulsive force from the star. 
























































^ we will discuss on the positions of the equilibrium points compared with the size of th4 

extended planetary atmospheres and show how the radiation pressure will influence the stability of the 
atmospheres of Hot Jupiters, with using HD 209458b as a case study. 


3. Results 

In this section [3l we present our results on the topology (13.1 p and on the decreasing size of the Roche 
lobe due to the radiation pressure. We derive an analytical formula to approximate the Hill’s sphere size 
depending on the radiation pressure intensity fl3.2p . We will then discuss the implications of the modihed 
Hill’s sphere on the planetary atmosphere stability fl3.3p . This modihcation will signihcantly affect the 
giant planets close to their host star such as the Hot Jupiters. This section will provide quantitative results 
for specihc situations including the case of HD 209458b (and for hydrogen atoms), but both the approach 
and the results can be essentially applied to any planetary atmosphere. 


3.1. Topology of equipotentials 

In this hrst part, we study the topology of equipotentials (or ZVC curves) modihed by the radiation 
pressure with the example of HD 209458b. 

The hgure [1] shows the topology of the equipotentials for increasing f3 values in the case of HD 209458b 
(i.e. /i ~ 5.83, see www.exoplanets.eu/HD209458b). Extremely close to its star (0.047 A.U.), with a mass 
(0.69 Mj) a nd a size (1.35 Rj) of the same o rder as Jupiter’s ones, HD 209458b is strongly affected by 


tidal forces flLecavelier des Etangs et al. 


20041) . Nevertheless, the photogravitational Three Body Problem 


was not yet covered concerning its implication on the stability of planetary exospheres; only the potential 
effects of tidal forces were investigated. We considered here a range of 0 values of 0 to 5, to be compared 


with t he ranges of 0 values for Hydrogen at HD 209458b as reported bv iBourrier and Lecavelier des Etansrs 


( 20131) ( these authors take into account the Doppler shift so that 0 depends on the velocity along the 
radial direction with respect to the star, /J(n 0) 4). 

As can seen in the hgure [H the different Lagrange points migrate in the direction toward the star. 
Indeed, the radiation pressure is a force opposed to the stellar gravity. For increasing 0 values, the 
equilibrium points need to be nearer to the star in order to put up the radiation pressure acceleration. 

For the critical value /J = 1, all Lagrangian points disappear except L 2 which should stay behind the 
planet. However, for any /J, L 2 becomes closer to the planet, and thus the sphere of inhuence or the Roche 
lobe has a decreasing size as a function of 0. 

In order to estimate the inhuence of 0 on the Roche lobe size, we derive in the next section a new Hill’s 


sphere radius taking into account the radiation pressure. 
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13 = 0 




(3 = 0.5 



(3 = 0.1 




(3 = 0.9 



Figure 1: Levels of effective potential for /3 = 0 (pure CR3BP, upper panels), 0.1, 0.5, 0.9, 1, 5 for HD 209458b. The 
crosses define the Lagrange points: Li (orange), L 2 (red), L 3 (blue), L 4 and L 5 (green). The central point corresponds to 
the star position. The black point around (1, 0) corresponds to the true planetary size and the surrounding white disk to the 
atmosphere with its supposed size. 
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We propose here to investigate the influence of the radiation pressure on the Roche lobe radius. As 
previously defined, the Roche lobe, in the case of = 0 (no radiation pressure), is the last equipotential 
closed around the planet and passing through Li. On another side, the Hill’s sphere radius is a first order 
approximation of the Lagrange Li and L 2 distances (Li and L 2 are approximatively at the same distance 

ill’s spher e. The radius of the Hill’s sphere is 


from the planet). Literally, the Roche lobe is inside the L 
approximately given at the first order by the value ^ /i/3 f Hill Il878h scaled by the star-planet distance. 


Following the demonstration for the Hill’s sphere radius Rh-, if we take into account the radiation 
pressure (i.e. /9 7 ^ 0), then Ru - i.e the position of the L 2 Lagrange point - is given by the following 
relation: 


-(1-/3) 


1 -p 


{1 + RhY 


_h_ 

Rl 


+ (1 — h T Rh) — 0 


—(1 — /3)(1 — fi)R^ — /i(l + Rh)'^ + Rh{^ + Rh)^{^ ~ h + Rh) — 0 


R^ T (3 — /RjR'ff -|- (3 — 2jR)R^ -f- [/3(1 — l^) ~ ^A^H — 2^R}{ — p — 0 


( 2 ) 


as presented in 


Simmons et al. 


fll985h fcf . equa tion 10 with 5\ = 1 — (3 and = 1 in their paper). 


According to Descartes’ rule (jPescartesl. 1163711 . this polynomial has only one positive root (because 


the sign between two successive monomials changes once), the Hill’s sphere radius. As we know, this 
polynomial does not have explicit roots (for polynomials with a higher de gree than 4, the roots cannot be 


written explicitly as a function of coefficients, cf. Abel-Ruffini theorem fj Abell . 


18241 ) L We can however 


reasonably assume that Rh 1, p 1 and p Rh- according to Hill’s formula (without radiation 
pressure), Rh{P = 0) = and according to flBeth et ahl. l2015bl) . a large radiation pressure leads 

to an exopause or Hill sphere radius located at Rh{/3) ~ (/i//3)^/^. We thus cut this relation at the first 
power not depending on f3 which is R?u (for the initial derivation of the Hill’s sphere radius, only the terms 
R?^ and 0(1) are kept). Thus, the relation [3^ gives us: 


?>R% + jSR^fj — fi ^ 0 


( 3 ) 


This polynomial has one positive root and two negative (if f3 > 9-^/i/12) or complex conjugates (if 
f3 < 9-^/i/12) roots. The value of the generalized Hill’s sphere radius is the positive root and can be found 
analytically thanks to the Cardano method: 

/3 ( 1 _, f 1 f9RHoV 


Rh = ^ \ 2 cosh | - argcosh y 2 \ j3 J 


( 4 ) 
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for f5 > 3-^9yu/4 with R^q = ^/JtJS the well known Hill’s sphere radius without radiation pressure. 


( 5 ) 


The upper panel of hgure [2] shows the evolution of the Hill’s sphere size as a function of /i and 13. The 
planetary systems with low mass planet/star mass ratios are thus more sensitive to the radiation pressure 
(in particular at low beta values , i.e. 0 to 0.5), with a Hill’s sphere radius decreasing more rapidly when 
the radiation pressure {(3) increases. 

To justify the use of this convenient formula, we compare our analytical solution of the equation |3] with 
the solution obtained numerically from the equation [2] (see hgure [2l middle panel). Our analytical formula 
thus approximates the position of L 2 with only ~ 3% error for jj, = 10“^ (case of the Sun-Jupiter system). 
Moreover, this error decreases quickly with /i decreasing. 

The Taylor series of our analytical formula gives at /3 = 0 ^ p/3 and at [3 —> -foo v^/V^. The 
13 = 0 asymptote corresponds to the classical Hill’s sphere radius without radiation pressure, whereas 
the {3 —> -|-oo asym ptote corresponds t o the exopause distance induced by the radiation pressure and 


discussed in detail by 


Beth et ah 


(l2ni5bh . Our generalized Hill’s sphere radius thus provides the location 


of the Hill’s sphere radius for any regime of the stellar radiation pressure. Mor eover, t h is giv es a math 


ematical basis for the radiat ion pressur e induc ed exopause hrst introduced by 


numerically demonstrated by 


Beth et al, 


BishopI fjl99lh and later 


fj2ni5bl) . The generalized formula is all the more interesting since 


the asymptotic values can be used for restricted conditions of the radiation pressure. For an error below 
~ 1% the hrst asymptote (i.e. the classical Hill’s sphere radius) is valid until (3 = O.IRho only, whereas 
the second asymptote (i.e. the exopause distance) is valid only above {3 = SORho- 


Planets 

/3 = 0 

(3 = 0.6 

/3 = 1.2 

/3 = 2 

/3 = 4 

/3 = 40 

Venus 

167.0 

36.22 

25.59 

19.80 

13.99 

4.42 

Earth 

234.9 

52.68 

37.23 

28.82 

20.37 

6.44 

Mars 

319.7 

49.63 

34.99 

27.06 

19.12 

6.04 

HD 209458b 

4.24 

2.14 

1.57 

1.23 

0.88 

0.28 


Table 1: Modified Hill’s sphere radius in planetary radii with the effect of the radiation pressure for different planets (Venus, 
Mars, Earth and HD 209458b) and for different /3 values (0,0.5,2,4,40). The classical Hill’s sphere radius (without radiation 
pressure) corresponds to /3 = 0. 

The f3 parameter depends both on the stellar photon hux at each wavelength and on the species 
considered. This work investigates only the radiation pressure effect on Hydrogen atoms but for future 
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Figure 2: (Upper panel) Ratio between the generalized Hill’s sphere radius and the planet-star distance as a function of /3. 
(Middle panel) Relative difference between the exact numerical root from the equation [2] and our analytical solution from the 
equation [3] (Lower panel) Same as upper panel, with the modified Hill’s sphere radius scaled with the Hill’s sphere radius 
for /3 = 0 (pure CR3BP), for different /r values. For /3 = 0, we have the usual value of the Hill’s sphere radius without 

radiation pressure. For higher 0 , the m odified Hill’s sphere radius converges asymptotically to y/ /i//3, i.e. the exopause 


distance discussed by 


Beth et al 


iaoisbl). 

















^or^"*?)ffi3'‘^pecies such as Helium should be investigated too. In the table [H we provide the modified 
Roche lobe size in planetary radii for Venus, Earth, Mars and HD 209458b for various ty pical d values. 


In the current Solar System, for Hydrogen, f3 is between 0.6 and 1.2 with daily variations (1 Vidal-Madiar 


f)l975l) ). The first column (i.e. beta=0) references the Hill’s sphere radius in planetary radii as used in 
literature. With the effect of the radiation pressure, the size of the Hill’s sphere may decrease strongly 
depending on the mass, the distance to the star and the radiation pressure (/5). For example, the size 
of the Hill’s sphere decreases by a factor of 6.5/6.3/9.1 for Venus/Earth/Mars when (3 = 1.2 (maximum 


radiation pressure) is considerec 
quite interesting: according to 


instead of d = 0 (no radiation pressure). The case of HD 209458b is also 


Bourrier and Lecavelier des Etangsl (120131) . the f3 value for HD 209458b 


should be ~ 4, neglecting relativistic effects, and thus the Hill’s sphere radius should be below the surface 
according to our study. 

We can also investigate the conditions encountered by the early Sola r System. The Lyman -g intensity 
was 20 times higher for the one million year old Sun than at present (iLammer et al.l. l2013l) . leading to 
an estimated (3 for Hydrogen of about 40. If we assume the planets were in the same place as today, the 
sphere of influence should be much closer to the planetary surface, at only a few planetary radii. 

The consequences on the atmospheric stability and escape flux, in particular for the case study HD 
209458b, will be more detailed in the next sections |33] and [TH 

3.3. Consequences on the stability of planetary atmospheres 


Regarding the example of the hot Jupiter HD 209458b, 


Lecavelier des Etangs et al 


(12004]) proposed 


that the exobase of the planet is shaped by the Roche lobe proximity, since the Li Lagrange point is 
expected to be close to the planet due to the close host star. This then induces a geometrical blow-off: 
the exospheric particles need much less kinetic/thermal energy to escape by reaching the close Roche lobe. 
Based on our analysis, we agree in this case with the geometrical blow-off of HD 209458b proposed by 


Lecavelier des Etangs et al 


(1200411 but for a different reason: as can be seen in figure |3l the radiation 


pressure pushes the position of the Roche lobe (i.e. Lf) much closer to the star and thus much further 
from the planet, but the L 2 point is pushed closer to the planet below the expected exobase location (and 
even below the “surface”), so that , no kinetic/thermal energy is needed for exospheric particles to escape 
even for a low exospheric temperature. 

As mentioned above the e arly Solar Syste i n plan etary atmospheres encountered a Lyman-a flux about 
20 times higher than today (ILammer et al.l. l2013l) and leading to a /J value of 40 instead of 0.6 — 1.2. 
Consequently, assuming the distance of Venus, Earth and Mars were the same as today, the sphere of 
influence for Hydrogen should be of the order of 4 and 6 planetary radii. The radiation pressure thus 
strongly limited the possible extension of the atmosphere: above this limit, the Hydrogen (and species 
with similar (3) could not be bounded to the planet by the gravity and easily escaped to the interplanetary 
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Figure 3: Equipotentials around HD 209458b for different radiation pressure intensities: /? = 0 (blue), 1 (green) and 4 (red). 
The corresponding Hill’s sphere is the circle in dotted line. Th e filled centered circle corre s ponds to the solid part of the 

(booj). 


planet and the surrounding circle to the exobase as proposed by 


Lecavelier des Etangs et al. 














iiectrttttiy-^-tiis radiation pressure effect is a new constraint for understanding the planetary atmospherfei 
evolution in the Solar or exoplanetary systems as well as for their modelling: the size of the simulation 
box should thus be at least larger than the Hill’s sphere radius of each species to appropriately take into 
account the radiation pressure effect. 

In the aim to have a complete picture of the atmospheric evolution of planets, which includes the 
influence of the radiation pressure, one should not only take into account the evolution of the stellar EUV- 
XUV fluxes that drive the heating/expansion of the atmospheres, but also the part of the flux spectrum 
that drives the radiation pressure. Depending on the flux spectrum of the host star, the species inside 


the atmosphere will b e more or less affected by t 


le ra diation pressure based on their specihc (3 value. For 


example, if we refer to lChamberlain and HnntenI fjl987l) for the Solar System conditions, the f3 parameter of 


Helium can be between 0 or very small values (e.g. 0.06 for He H) and 260 depending on the wavelength. If 
several species are not equally sensitive to the stellar radiation pressure, this could then affect the dynamic 
of each species above the exobase and yield to a differential drag of the species in the upper atmosphere, 
or to a protection of (radiation pressure) sensitive species such as Hydrogen by less sensitive ones through 
collisions. This last process could help understanding how Hydrogen rich atmospheres can remain stable 
for long periods despite a blow-off type escape regime expected. These questions will be investigated in a 
future work. 

Moreover, it can be reasonably assumed that the radiation pressure will have a different effect on 
planetary atmospheres according to the Knudsen number Kn\ 


• for low Kn (collisional regime), i.e. the scale height is greater than the mean free path, the radiation 
pressure will increase the mechanical energy of each Hydrogen atom on short distances similar to the 
mean free path. Then, through collisions, this increase will be redistributed quickly to other species 
and transformed into heat. The escaping flux will be a function of the incoming UV flux and its 
heating efficiency. 


• for high Kn (non collisional regime), i.e. the mean free path is higher than the scale height, the 
radiation pressure will affect only the dynamics of each individual species on time scale shorter than 
1/z/, z/ the collisional frequency. The escaping flux will be function of the radiation pressure and the 
local thermal speed. 


3.4- Escaping flux: the case study HD 209458b 

As detailed above (and seen in table [T]), the hot Jupiter HD 209458b has its Hill’s sphere radius 
(or exopause) below the exobase and even below the surface due to the strong radiation pressure of the 
host star. Since the exonanse is below the exobase, all particles at the exobase thus escape fsee Beth 

et al. i2015b H. This means the escape velocity is “virtually” reduced to 0 and thus the Jeans’ parameter 













^esc/ ^th K^esc aiLyjL Vth being the escape and ther mal speed) ten ds to 0. The critical Jeans’ paramet^ 
for the hydrodynamical blow-off regime proposed by iHnntenI (119821 1 is 2, but with the additional effect of 
the radiation pressure, Ac will reach even lower values. The thermal escape flux may then be given by the 
Jeans’ escape formula applied to the whole exobase with Ac = 0, which leads for Hydrogen at HD 209458b 
to the following mass loss rate: 


M ^ AT:rl^jnFj{\c ^ 0) ^ ~ 5.10^‘’g.s ^ 


( 6 ) 


for a temperature T^xo of 8000 K , an exobase at a distance o 
exoba se of the order of 10^^ m“^ flLecavelier des Etangs et ah 


2.8 p: 


(2004h 


(2013)). This loss rate is in agreement with the estimates by 


anetary radii and a density at the 


Bourrier and Lecavelier des Etangs 


Vidal-Madiar and Lecavelier des Etansrs 


(120041 1. i.e. above 10® g.s and of the sarne orde r as recent models which predict a range of loss rates 


and could be much lower than the one used by 


between 10® and 10^^ g.s ^ (IKoskinen et al.l ( 2013111. However, the exobase tem perature is poorly known 


Lecavelier des Etangs et al 


(1200411 to match their model with 


the spectroscopic observations. Then, we have perform ed as well the same c alcula t ion with the effecti ve 

( 12000 ): 


temperature of HD 209458b which is T^ff ~ 1400 K (iCharbonneau et al 


Evans et al. 


(1201,511 1 


Assuming the same order of magnitude for the density at the exobase than previously assumed, we obtain 
M = 2.10^® to be compared with the Jeans’ escape, which is here of about M 20 g.s-h 

We thus obtain similar escape rates to previous models, but with a different process, i.e. with a 
geometrical “blow-off’ due to the radiation pressure influence and not due to the Roche lobe proximity. 
Furthermore, for Hot Jupiters, the thermal escape or Jeans’ escape is in the literature neglected but our 
study shows that the radiation pressure acceleration modihes the gravitational potential in such a way 
that the exospheric particles are not bounded any more and/or the Jeans’ parameter tends to 0. Thus, 
the thermal escape is strongly enhanced by the stellar radiation pressure. Finally, our study shows that 
extreme temperatures, compared with Te//, are not required to be in agreement with the expected mass 
loss rates for HD 209458b. 


4. Conclusions 

In this paper, we show the effect of the radiation pressure on the Circular Restricted Three Body 
Problem (CR3BP). This problem was previously tackled for Celestial Mechanics purposes but rarely in 
the context of Planetary Science. We hrst discuss the influence of the radiation pressure on the positions 
of the Lagrange points, before we derive for the hrst time an analytical formula to approximate the Hill’s 
sphere radius with the additional effect of the radiation pressure. We highlight the strong interest of this 
derivation for the case studies of exoplanets and in particular Hot Jupiters. Most of the papers dealing 
with the Hill’s sphere (or Roche lobe) and the stability of Hot Jupiters atmospheres do not take account 















































pressure on the CR3BP. We show in extreme cases such as for Hot Jupitel:^ 
like HD 209458b that the Hill’s sphere radius could be drastically reduced and even be located below 
the surface. The components of the upper atmosphere, which are sensitive to the radiation pressure, are 
thus not “bounded” to the planet any more and then escape to the interplanetary medium. The resulting 
increased escape calculated for HD 209458b is in agreement with the observations. Moreover, we show 
that this Hot Jupiter is in a “blow-off’ type escape regime due to the radiation pressure of the host star, 
in contradiction with previous works that considered a blow-off regime due to the proximity of the Roche 
lobe (the Roche lobe being actually pushed toward the star by the radiation pressure). Thus, the radiation 
pressure can modify the planetary escape flux and strongly affect the atmospheric evolution of hot Jupiters 
atmospheres from their early age, depending on the evolution of the host star as well. 

In a future work, we will investigate how the intense radiation pressure of the young Sun may have 
influenced the atmospheric escape and evolution of the inner Solar System planets, whose Hill’s sphere 
radius was significantly reduced. 
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